Acoustic oscillations and viscosity 



Roy Maartens* 

School of Computer Science and Mathematics, Portsmouth University, Portsmouth P01 2EG, Britain 

Josep Triginert 

International University of Catalonia, C/Immaculada 22, 08017 Barcelona, Spain 
and Departament de Fisica, Universitat Autonoma de Barcelona, 08193 Bellaterra, Spain 

(July 1998) 



Using a simple thermo-hydrodynamic model that respects relativistic causality, we revisit the 
analysis of qualitative features of acoustic oscillations in the photon-baryon fluid. The growing pho- 
ton mean free path introduces transient effects that can be modelled by the causal generalization 
of relativistic Navier-Stokes-Fourier theory. Causal thermodynamics provides a more satisfactory 
hydrodynamic approximation to kinetic theory than the quasi-stationary (and non-causal) approx- 
imations arising from standard thermodynamics or from expanding the photon distribution to first 
order in the Thomson scattering time. The causal approach introduces small corrections to the 
dispersion relation obtained in quasi-stationary treatments. A dissipative contribution to the speed 
of sound slightly increases the frequency of the oscillations. The diffusion damping scale is slightly 
increased by the causal corrections. Thus quasi-stationary approximations tend to over-estimate 
the spacing and under-estimate the damping of acoustic peaks. In our simple model, the fractional 
corrections at decoupling are > 10 -3 . 
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Acoustic oscillations in the photon-baryon fluid before decoupling leave a vital imprint on the cosmic microwave 
background [jjj (and possibly also on large-scale structure The form of this imprint encodes information on fun- 
damental cosmological parameters, and increasingly accurate and refined numerical integrations are being performed 
to produce predictions that can be tested against current and upcoming observations. As a complement to detailed 
numerical models, it is also useful to analyze qualitatively and analytically the key physical features such as acoustic 
oscillations. As pointed out by Hu and Sugiyama [Q, numerical integrations are sufficient for direct comparison of 
specific models with observations, but do not readily produce a qualitative and analytic understanding of the physical 
processes at play. In ffl , acoustic oscillations and their damping due to photon diffusion are analyzed analytically via 
^ ■ expanding the integrated Boltzmann multipoles in the Thomson scattering time r T . To zero order in r T , i.e. in the 
tight-coupling approximation, which is valid on scales much larger than the photon mean free path, the oscillations 
are undamped. Damping arises from the first order approximation, which introduces a shear viscosity via the radi- 
ation quadrupole. This approximation is hydrodynamic, in the sense that all multipoles beyond the quadrupole are 
neglected, and it is effectively equivalent (see S) to the non-equilibrium thermodynamic approach of Weinberg [Q, 
which is based on the relativistic Navier-Stokes-Fourier theory developed by Eckart. 

As pointed out by Peebles and Yu expansions in t t , even beyond first order (see 0), cannot take proper 
account of the change in the photon mean free time, especially near to decoupling. In other words, the 0(t t ) 
approximation, and the equivalent Navier-Stokes-Fourier approximation, are inherently quasi-stationary, and provide 
a limited approximation to the Boltzmann equation. An improved hydrodynamic approximation to kinetic theory is 
the causal thermodynamics developed by Israel and Stewart which generalizes the non-causal Eckart theory by 
incorporating transient non-quasi-stationary effects. In this paper, we use causal thermodynamics to refine aspects 
of the work by Weinberg Q and Hu and Sugiyama Of course, all hydrodynamic approximations, which assume 
that the photon-baryon fluid is close to equilibrium via interactions, will break down when the Thomson interaction 
rate becomes too low, i.e. as decoupling is approached. 
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In a previous paper |J7j , we developed a general formalism for incorporating causal thermodynamics into perturba- 
tion theory, thus providing a self-consistent approach to density perturbations in dissipative cosmological fluids. Our 
formalism is based on the covariant perturbation theory of Ellis and Bruni ||||, and the covariant causal thermo- 
dynamics of Israel and Stewart [ra . Here we apply the general formalism to the case of viscous damping of density 
perturbations in the photon-baryon fluid. This effect has been comprehensively analyzed via detailed study of the 
Boltzmann equation (see e.g. |5|Jl0|-p^] ) . Dissipative hydrodynamics provides a simplified alternative to a full kinetic 
theory analysis, that can illuminate some of the key physical effects without the same level of detail and complexity. 
The approaches of Weinberg Q and Hu and Sugiyama [Q neglect gravitational effects, which is not unreasonable on 
subhorizon scales. We refine their approach by including metric perturbations. More importantly, their approach is 
inherently quasi-stationary, and the main refinement we introduce is to incorporate transient effects via the causal 
generalization of Navier-Stokes-Fourier thermodynamics. 

It turns out that the relaxational effects which are incorporated into the causal transport equation, but which are 
neglected in the non-causal theory, have a small but interesting impact. The sound speed, which determines the 
frequency of acoustic oscillations for each mode, acquires a positive dissipative correction. This produces a small 
increase in the frequency of acoustic oscillations on each scale. The viscous damping rate found in quasi-stationary 
approximations [EM] also acquires a positive correction, leading to a small increase in the diffusion cut-off scale. A 
rough estimate based on our simplified model is that the fractional corrections are > I0 -3 by decoupling. 

For simplicity, since we are interested in qualitative features rather than specific numerical predictions, and since 
we wish to focus on relaxational corrections to quasi-stationary models, we assume a flat background universe with 
negligible non-baryonic content. We treat the mixture of baryonic matter and radiation as a single dissipative fluid, 
with dissipation arising from the growing mean free path of photons in their Thomson scattering off electrons. We 
consider the era between matter-radiation equality and decoupling, neglecting the matter pressure. Dissipative effects 
are expected to be greatest near to decoupling. We will also neglect thermal conduction for simplicity, while the bulk 
viscous stress is necessarily negligible since the fluid is in the non-relativistic regime. Thus the dissipative effect of 
photon diffusion is reduced (within a radiative hydrodynamic model) to a shear viscous stress n a b- 

In non-causal thermodynamics (and in an 0{j T ) expansion of the the Boltzmann multipoles), this shear viscous 
stress is algebraically determined by the shear a a b via [Q 

TT ab = -2rj(Jab , (f ) 

where r\ is the viscosity. The implicit instantaneous relation between the cause (shear, i.e. anisotropic expansion rate) 
and the effect (anisotropic stress), is what underlies the pathologies of the theory, i.e. that it admits dissipative signals 
with infinite wavefront speeds and that its equilibrium states are generically unstable |l3| ]. In the causal theory, the 
stress is no longer algebraically and 'instantaneously' determined by the shear, but satisfies an evolution equation in 
which there is a time-lag between cause and effect ||: 

T-kab + K ab = -2T)O ab , (2) 

where r is the relaxation time-scale. This causal transport leads to a subluminal speed of viscous pulses, and introduces 
small but interesting corrections to the quasi-stationary dispersion relation. 

Causal thermodynamics has a solid kinetic theory motivation via the relativistic Grad I4-moment approximation 
to the Boltzmann equation This approximation effectively restores second-order non-equilibrium terms in the 
entropy that are neglected in the relativistic Chapman-Enskog approximation |fl4j| , which leads to the non-causal 
theory. Thus causal thermodynamics provides a more satisfactory and less incomplete hydrodynamic approximation 
to kinetic theory than the standard non-causal theory. The differences are greatest for high frequency phenomena, on 
scales comparable to the mean free time or path, when the Chapman-Enskog approximation breaks down, while the 
Grad approximation remains valid. 

We follow the notation of Q|. Units are such that c = 1 = 8ttG, except where numerical values are given in specified 
units. The present Hubble rate is Hq = fOO/i km/s/Mpc. 



II. CAUSAL VISCOUS PERTURBATIONS 

Ref. [Q contains a general discussion and general covariant perturbation equations for causal dissipative hydrody- 
namics in cosmology. In the interests of brevity, we will not repeat here the derivations and motivations of [Q , but 
simply quote the results as they are needed. 

In the era between matter-radiation equality and decoupling, the mixture of photons, baryons and electrons is 
treated as a single dissipative fluid, with dissipation due to photon Silk diffusion Jl5|]. The baryons and electrons are 
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tightly coupled, culminating in their recombination, while the photon coupling to baryonic matter through Thomson 
scattering weakens as the mean free path grows. We assume for simplicity that there is no non-baryonic cold dark 
matter. We also assume negligible thermal conduction and particle flux , so that the particle and energy frames 
coincide Q, and we can choose the fluid four-velocity u a (where u a u a = — 1) as the four-velocity of this frame. The 
bulk viscous stress is negligible, since the fluid is non-relativistic. 

The fundamental quantity in the covariant approach to density perturbations H is 



aD a 6 a 



,D 2 p 



(3) 



where D a is the covariant spatial derivative, a is the background scale factor, and 



is the comoving fractional energy density gradient, which is a covariant measure of density inhomogeneity. The total 
energy density p — p B + p a is made up of baryonic and radiation parts, and 5 is the density fluctuation scalar for the 
photon-baryon system, considered as a combined single fluid. 

Covariant entropy perturbations in a single dissipative fluid model (as opposed to a 2-component model), are 
sourced by heat flux and bulk viscous stress, and shear viscous stress has no effect to linear order (Eq. (18), 0). If, 
as we assume, the initial entropy perturbation is zero, there are therefore no entropy perturbations.n Then S satisfies 
the evolution equation (Eq. (28), 0) 



5 + 2H(l- 3w + §c 2 ) 6 - \H 2 (1 + 8w - 3w 2 - 6c 2 ) 5 - c 2 s D 2 5 = S[tt] , 



(4) 



where an overdot denotes the covariant derivative along u a , H = a/a, w = p/p, with p = p B + p R the total isotropic 
pressure, and c 2 = p/p is the adiabatic sound speed. The shear viscous source term is given by (Eq. (32), 0) 



S[tt] = 3HS - 3iJ 2 (l + 6w - 3c 2 )<S + D 2 <S , 



(5) 



where the covariant shear viscous stress scalar is 



S = a 



2 D a D b 7r Q6 



The viscous transport equation 
tS 



leads to the following evolution equation for S (Eq. (46), 0):0 



l-H I 3(1 + w)t 



4?7 



4?/ 



3p(l + w) 



5 - 3wH5 



(6) 



Equations (f§)-(||) form a coupled system that governs the evolution of density perturbations with causal viscosity. 

The thermal energy of matter is much less than the rest mass energy after the epoch t e of matter-radiation equality, 
so that p B ~ 0, and then p » \p R . Taking a = 1 at the present epoch, we have in the background (where matter and 
radiation are non-interacting) 
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W= 3 
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21 ' 



(7) 



for a > a c . To a reasonable approximation (sufficient for the purposes of our simple model), we can neglect these 
quantities relative to 1 in the round brackets in equations (|])-(|^). Furthermore, since w -C 1 and rjH/ p <C 1, we can 
neglect the last term on the right of Eq. (^) . 

We decompose into Fourier modes 8(t, x) — > A(t, k) and S(t, x) — > £(t, k), where k is the comoving wave number, 
so that the proper wave number is K = k/a. The proper wavelength A = 2t:/K is constrained by A T < A < X H , 



In any case, the entropy perturbations are decoupled from the density perturbations, because the non-barotropic index 
(dp/ds) p vanishes in the background, so that the entropy source term in the density perturbation equation is zero (Eq. (28), 



i) 



The p factors are mistakenly omitted in equations (46) and (54) of 
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where X H = H 1 is the Hubble scale (above which thermo-hydrodynamic effects do not operate), and the minimum 
scale X T is the photon mean free path for Thomson scattering: 



A^ 



1 



(8) 



Here a T is the Thomson cross section and n is the free electron number density. Perturbations on scales below 
A T will be wiped out by photon streaming (and the hydrodynamic model breaks down on these scales). Causal 
thermodynamics applies down to the Thomson scale, but the non-causal theory requires A ^> A T . The damping effect 
of photon streaming operates beyond the Thomson mean free scale, since photon diffusion out of over-dense regions 
rapidly destroys perturbations below a critical scale, as shown by Silk fl5|| . The critical cut-off scale is the diffusion 
scale A D , which is considerably greater than the Thomson mean free scale, as we will confirm below. 

Collecting the above points, and using the background Friedmann equation p — 3-ff 2 , the coupled system becomes 



A • 



2H A - \H 2 (1 - § c 2 K 2 ) A = 3H± - 3H 2 (l -I- \K 2 ) S , 



3Ht± + 3H(l 



77* — 3Ht) £ = ?7*A . 



where we have defined the dimcnsionless expansion-normalized shear viscosity 

-ir] 



V* = 



3H ' 



and proper wave number 



The relevant scales constrain K* by 



K K 0^ 



(9) 
(10) 



(11) 



(12) 



2vr < K* < 2iT-Z- 



(13) 



Radiative shear viscosity due to Thomson scattering is given by (see and ||16|-|18| for refinements using kinetic 
theory 3 ) 



(14) 



where ro is the blackbody radiation constant, T is the photon temperature and r T = A T is the photon mean free time 
for Thomson scattering. The causal relaxation time-scale r is by Eq. (||) the characteristic time taken for the fluid to 
relax toward equilibrium if the viscous 'driving force' were to be 'switched off'. Thus we expect r is of the order of a 
few Thomson interaction times, and below the characteristic diffusion time: 



< 



T < T r 



(15) 



where t d = A D . 



In the era we are considering, H decays approximately like a 3 / 2 , n decays faster than a 3 because of recombination 
effects, and T decays as a -1 . Thus to lowest order 



K* = (#*) e 



1/2 



3/2 



(16) 
(17) 



3 Note that claims a correction factor of ^ in Weinberg's expression for 77. This arises from anisotropic scattering effects, 
and a further small correction is also induced by polarization effects jjj. These corrections are not important for our simple 
model. 
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3+1 



3/2+1 



1/2 



(18) 
(19) 
(20) 



where 7 > encodes the growing effects of recombination on the number density of free electrons, the details of 
which are not important for our simple model. Using the standard numerical values Jis|| of t c rj 3 x 10 10 / 
redshift z c «2x 10 4 /i 2 , and present baryon number density n E 
luminescent) density fraction f2 B « 1, we find that 



10 8 h 8 cm 3 , with the present baryonic (dark and 



1.5 x 10 21 /i" 4 cm, (A, 



10 ib /T 8 cm, (n* 



1.5 x 10~ 6 h~ 4 . 



(21) 



Equations (^|) and ( jio| ) govern the evolution of density perturbation modes within a simplified causal viscous fluid 
model. The system may be decoupled to produce a third-order equation in A. In practice we analyse the coupled 
system, but for completeness, the decoupled equation is 

r A + [1 + \tH + t + F] A 
+ 77 [2 + (| + ±7f 2 ) jj, - f r77 + 2f - (tH)- 1 {rr?* + F(t?* - 2r77)}] A 



77^ 



where 



(1 + 77* - \tH + r + 7 1 ) (1 - \c\Kl) -tH- \rKl (c s 2 ) ' 



F = —t [In |r77 (6 - K f) - 3(1 + ^)|] ' , 



A = 0, 



(22) 



and we used the background field equation 77 = — 177 2 . In the non-causal case r = 0, Eq. 



3) reduces to the second 



order equation 



A + 77 



V* 

77(l+?7* 



A-§77 2 [(l + ^)(l-fc 2 K 2 )]A = 



(23) 



III. VISCOUS ACOUSTIC OSCILLATIONS 



The key features of the solutions are readily analyzed qualitatively, given the simplicity of our model. We assume 
that the coefficients in equations @ and ( |i"o| ) are slowly varying, which is reasonable for small-scale modes. Then 
we can try a WKB solution of form A(t,k) = a~ 1 A(k) expi fuidt and k) = a^ 1 B(k) expi J udt, where the 
scale factor terms remove the non-dissipative effects of expansion on small-scale modes, facilitating comparison with 
Weinberg's results, which neglect expansion. This gives 





'A' 




B 



0. 



±77 2 (5 - IclKl) -ZiHu + (6 + Kf)H 2 

777*0; + i]*H SirHuj + 3(1 + 77* — 4tH)H 

Non-trivial solutions require a vanishing determinant of the coefficient matrix, leading to the dispersion equation 

2iruj 3 + 2(1 - 4t7/> 2 - i [277* ( 2 + \ K *) + rH ( 2c I K * - 5)] Huj 
+ [(1 + r;* - 4r77) (5 - 2c 2 s Kl) + 2^ (l + ±7^ 2 )] 77 2 = . (24) 

In general, the solutions will have the form 

uo = ±Kv + i T , (25) 

where v is the (real) phase speed of oscillations and T is the viscous damping rate. Damping will be severe for scales 

A < A D , where the diffusion cut-off scale is defined by equality of the damping rate and the expansion rate, i.e. 
r(AJ=77. 
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In the dispersion relation ( pif ) , we can neglect 77* relative to 1 by equations (^) and (^Tj) . Equations (Q) , (|2|) , ( |T7| ) 
and (|T]) show that for the relevant small scales (A <C A H ), we havey Cgif: 2 ^> 1. Then the dispersion cubic may be 
written 



it^I + (1 - At^ujI - 1 (±77*X 2 + T.clKl) «,-(!- 4t*) c s 2 ]f, 2 = . 



(26) 



where = oj/i/ and t* = t/H 1 are dimensionless expansion-normalized variables. For the non-causal theory, with 
t = 0, Eq. (pq) becomes 



3w»- (i^)^-3cX 2 = 0, 
where an overbar indicates a quantity in the non-causal theory. Eq. (^) has solution 

f 2 



(27) 



-2 2 



1 



2 K 2 



Since T/H < 1 for the scales that survive viscous damping, and c 2 K 2 3> 1 for the relevant sub-horizon scales, we 
have v 2 = c 2 to a good approximation. Thus 



r = k- 



(28) 



which agrees (when the isotropic pressure, heat flux and bulk stress are negligible) with Weinberg's special relativistic 
generalization of the non-relativistic Landau-Lifshitz result B .H The diffusion cut-off scale follows from Eq. (ES) as 



5 - 2lT 211 x 



(29) 



By equations (|f7j), ( p!8| ) and (|20|), the diffusion scale is well above the Thomson scale, although the ratio decreases 
with expansion: 



A,, V a / 



5/4+/ 



(30) 



Now we consider the causal corrections to Eq. (p8|). From equations (|l5|), (|l9|), ( pf| ) and (p9|), it follows that 
r* = <C 1. We write = tD* (1 + yr*), where satisfies Eq. (p7j), and solve the dispersion cubic Eq. ( |26| ) to first 
order in r*. This gives 



H 



(f + AH) 
±c s K 



on using c 2 if 2 > 1. It follows that 



v = c s (1 + rf ) , 

r = f [1 + 2t (f + 2H)] 



(31) 
(32) 



Thus the causal generalization of standard thermodynamics leads to corrections that increase the effective sound 
speed and damping rate. The latter increase leads to an increase in the diffusion cut-off scale, which is determined 
by setting f » H in Eq. <M): 



X D (I + 3tH) 



(33) 



In order to estimate the size of the corrections, we need an expression for the causal relaxation timescale t. In 
the case of a simple (i.e. single-species) relativistic gas (obeying classical or quantum statistics), the thermodynamic 
coefficients such as r and rj are in principle known via complicated kinetic theory formulas involving integrals over 



Numerical integrations in the special case c\Kl < 1 are given in Q . 
5 See also [^l) for a kinetic theory approach, and |22j for the inclusion of magnetic fields. 
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energy and collision cross-sections ||. In the radiative transfer context, where 2 particle species are involved, one 
treated hydrodynamically, the shear viscosity is given by Eq. (jlj), while r is implicitly assumed to be zero in quasi- 
stationary treatments p|, [i^ , |l7| . If Grad-type moment approximations are used, then r is equal to the mean collision 
time as a consequence of the hydrodynamic assumptions |^,^|. However, it seems physically reasonable that r could 
be greater than t t , if the hydrodynamic assumptions in radiative transfer are somewhat relaxed. Thus a simple 
generalization 

r = olt t , (34) 

where a > 1 is a constant, seems reasonable within the phenomenology of thermodynamics. 

The correction to the diffusion scale, and the maximal correction to the sound speed (which occurs at T = H), are 
determined by tH. By equations (0), ( |l9|) and (pl|), we find that in our simple model 



\ \ / \ 3/2+/ / x 3/2+7 



a. 



tH = a[^) [ — ) « 7a x lO^/i" 4 ( — ] . (35) 



The correction factor grows with expansion, as expected, since the mean free photon path is growing. We take h = 0.6. 
At matter-radiation equality, the fractional correction is then w5ax 10~ 5 >5x 1CP 5 . As decoupling is approached, 
the ionization parameter / becomes more important. However, for a rough under-estimate, we can set 1 = 0. (In 
any case, all hydrodynamically based models will break down near decoupling.) With z<j oc ss 1100, we find that the 
maximal fractional correction introduced by relaxational effects is 

(^)max > " x !0~ 3 £ 10~ 3 , (36) 
in our simple model. This is the relevant estimate, since the imprint of acoustic oscillations is frozen in at decoupling. 



IV. CONCLUSIONS 



Our principal result is the causal generalization and upward correction of the quasi-stationary sound speed, given 
in Eq. ([H]), and damping rate, given in Eq. (|32|). These corrections are small, at roughly the 0.1-1 % level, as 
given in Eq. (|3^), but they show how quasi-stationary approximations tend to under-estimate the frequency and 
damping of acoustic oscillations, and they indicate the refinements that would arise from a more complete kinetic 
theory approach, which involves far greater computational complexity. The qualitative features of the refinements 
are: (a) a higher threshold for the survival of density perturbations, i.e. an increased diffusion cut-off scale, given by 
Eq. ([53() ; (b) the appearance of a dissipative contribution to the sound speed of the photon-baryon fluid (consistent 
with the general theoretical results for perturbations of stationary fluids in flat spacetime |]|b|) which increases the 
frequency of oscillations; (c) the scale-dependence of the corrections to the damping rate and sound speed; (d) the 
central role played by the relaxation time-scale r, which provides a simple scalar encoding for the non-quasi-stationary 
effects of photon diffusion. 

This paper shows that in principle dissipative hydrodynamics with causal transport equations can provide an 
economical alternative to kinetic theory models, with some additional insights into the nature and source of acoustic 
oscillations. The simplified model used here could be refined, principally by incorporating thermal conduction, to 
produce an improved model that nevertheless will still allow a qualitative and analytic treatment, unlike the full 
kinetic theory approach. 

Our model can also be seen as part of a growingbody of examples in cosmology (e.g. p4j), relativistic astrophysics 
(e.g. P5| ), and other areas of physics (e.g. ff(|)£] which show that interesting and sometimes significant physical 
differences can arise from the causal approach to thermodynamics. (A general overview is given in p7|.) Even if 
the causal corrections are small, as is the case here, they provide further insight into the physics, with the added 
advantage that one avoids the essentially unsatisfactory features of non-causal thermodynamics. The predictions of 
the non-causal theory are readily recovered in the appropriate limit. 



In non-relativistic physics, causal theories are usually called 'extended' or 'hyperbolic'. 
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